Differential selection of yield and quality traits has shaped genomic signatures of cowpea domestication and improvement.

Cowpeas (tropical legumes) are important in ensuring food and nutritional security in developing countries, especially in sub-Saharan Africa. Herein, we report two high-quality genome assemblies of grain and vegetable cowpeas and we re-sequenced 344 accessions to characterize the genomic variations landscape. We identified 39 loci for ten important agronomic traits and more than 541 potential loci that underwent selection during cowpea domestication and improvement. In particular, the synchronous selections of the pod-shattering loci and their neighboring stress-relevant loci probably led to the enhancement of pod-shattering resistance and the compromise of stress resistance during the domestication from grain to vegetable cowpeas. Moreover, differential selections on multiple loci associated with pod length, grain number per pod, seed weight, pod and seed soluble sugars, and seed crude proteins shaped the yield and quality diversity in cowpeas. Our findings provide genomic insights into cowpea domestication and improvement footprints, enabling further genome-informed cultivar improvement of cowpeas.

Genome-wide association analysis reveals the optimal genomic regions for pod size in bean.

The snap bean is the most commonly grown vegetable legume worldwide, and its pod size is both an important yield and appearance quality trait. However, the improvement of pod size in snap beans grown in China has been largely hindered by a lack of information on the specific genes that determine pod size. In this study, we identified 88 snap bean accessions and evaluated their pod size traits. Through a genome-wide association study (GWAS), 57 single nucleotide polymorphisms (SNPs) significantly associated with pod size were detected. Candidate gene analysis showed that cytochrome P450 family genes, WRKY, and MYB transcription factors were the predominant candidate genes for pod development, and eight of these 26 candidate genes showed relatively higher expression patterns in flowers and young pods. A significant pod length (PL) SNP and a single pod weight (SPW) SNP were successfully converted into kompetitive allele-specific polymerase chain reaction (KASP) markers and validated in the panel. These results enhance our understanding of the genetic basis of pod size, and also provide genetic resources for the molecular breeding of pod size in snap beans.

Genome-wide characterization and expression analysis of the MLO gene family sheds light on powdery mildew resistance in Lagenaria siceraria.

MLO (mildew locus O) genes play a vital role in plant disease defense system, especially powdery mildew (PM). Lagenaria siceraria is a distinct Cucurbitaceae crop, and PM is one of the most serious diseases threatening crop production and quality. Although MLOs have been exploited in many Cucurbitaceae species, genome-wide mining of MLO gene family in bottle gourd has not been surveyed yet. Here we identified 16 MLO genes in our recently assembled L. siceraria genome. A total of 343 unique MLO protein sequences from 20 species were characterized and compared to deduce a generally high level of purifying selection and the occurrence of regions related to candidate susceptibility factors in the evolutional divergence. LsMLOs were clustered in six clades containing seven conserved transmembrane domains and 10 clade-specific motifs along with deletion and variation. Three genes (LsMLO3, LsMLO6, and LsMLO13) in clade V showed high sequence identity with orthologues involved in PM susceptibility. The expression pattern of LsMLOs was tissue-specific but not cultivar-specific. Furthermore, it was indicated by qRT-PCR and RNA-seq that LsMLO3 and LsMLO13 were highly upregulated in response to PM stress. Subsequent sequence analysis revealed the structural deletion of LsMLO13 and a single nonsynonymous substitution of LsMLO3 in the PM-resistant genotype. Taken all together, it is speculated that LsMLO13 is likely a major PM susceptibility factor. The results of this study provide new insights into MLO family genes in bottle gourd and find a potential candidate S gene for PM tolerance breeding.

Elucidation of the Flavor Aspects and Flavor-Associated Genomic Regions in Bottle Gourd (Lagenaria siceraria) by Metabolomic Analysis and QTL-seq.

Bottle gourd (Lagenaria siceraria) is a commercially important cucurbitaceous vegetable with health-promoting properties whose collections and cultivars differ considerably in their flavor aspects. However, the metabolomic profile related to flavor has not yet been elucidated. In the present study, a comprehensive metabolite analysis revealed the metabolite profile of the strong-flavor collection "J120" and weak-flavor collection "G32". The major differentially expressed metabolites included carboxylic acids, their derivatives, and organooxygen compounds, which are involved in amino acid biosynthesis and metabolism. QTL-seq was used to identify candidate genomic regions controlling flavor in a MAGIC population comprising 377 elite lines. Three significant genomic regions were identified, and candidate genes likely associated with flavor were screened. Our study provides useful information for understanding the metabolic causes of flavor variation among bottle gourd collections and cultivars. Furthermore, the identified candidate genomic regions may facilitate rational breeding programs to improve bottle gourd quality.

Unravelling the Genetic Architecture of Rust Resistance in the Common Bean (Phaseolus vulgaris L.) by Combining QTL-Seq and GWAS Analysis.

The common bean (Phaseolus vulgaris L.) is the most important legume crop directly used for human consumption worldwide. Bean rust, caused by Uromyces appendiculatus, is a devastating disease and usually causes severe loss of seed yield and pod quality. Deployment of resistant cultivars is the best strategy to combat this disease. However, despite being the largest snap bean-producing country, the genetic basis research of rust resistance has largely lagged in China. In this study, an RIL population and a diversity panel were evaluated for rust resistance against a purified rust isolate Cua-LS using a detached leaf assay. Deploying a QTL-Seq analysis in the RIL population, a 1.81 Mb interval on chromosome 4, a 2.73 Mb interval on chromosome 5 and a 1.26 Mb interval on chromosome 6 were identified as major QTLs for rust resistance, designated as Qur-1, Qur-2 and Qur-3, respectively. Through a GWAS diversity panel, 64 significant SNPs associated with rust resistance were detected, distributed in all 11 chromosomes and explaining 19-49% of the phenotypic variation. Synteny analysis showed that Qur-2 was validated in GWAS, but the rust QTL/SNPs detected in our study were different from the known genes, except Ur-11. A total of 114 candidate genes, including the typical NBS-LRR genes, protein kinase superfamily proteins and ABC transporter family proteins, were identified and proposed as the likely candidates. The identified 17 resistant accessions will enrich the resistant germplasm resources, and the detected QTLs/SNPs will facilitate the molecular breeding of rust resistance in the common bean.

Insight into the bZIP gene family in Lagenaria siceraria: Genome and transcriptome analysis to understand gene diversification in Cucurbitaceae and the roles of LsbZIP gene expression and function under cold stress.

The basic leucine zipper (bZIP) as a well-known transcription factor family, figures prominently in diverse biological and developmental processes and response to abiotic/biotic stresses. However, no knowledge of the bZIP family is available for the important edible Cucurbitaceae crop bottle gourd. Herein, we identified 65 putative LsbZIP genes and characterized their gene structure, phylogenetic and orthologous relationships, gene expression profiles in different tissues and cultivars, and responsive genes under cold stress. The phylogenetic tree of 16 released Cucurbitaceae plant genomes revealed the evolutionary convergence and divergence of bZIP family. Based on the specific domains, LsbZIP family were classified into 12 clades (A-K, S) with similar motifs and exon-intron distribution. 65 LsbZIP genes have undergone 19 segmental and two tandem duplication events with purifying selection. The expression profiling of LsbZIP genes showed tissue-specific but no cultivar-specific pattern. The cold stress-responsive candidate LsbZIP genes were analyzed and validated by RNA-Seq and RT-PCR, providing new insights of transcriptional regulation of bZIP family genes in bottle gourd and their potential functions in cold-tolerant variety breeding.

Identification of Bottle Gourd (Lagenaria siceraria) OVATE Family Genes and Functional Characterization of LsOVATE1.

The OVATE gene family is a class of conserved transcription factors that play significant roles in plant growth, development, and abiotic stress, and also affect fruit shape in vegetable crops. Bottle gourd (Lagenaria siceraria), commonly known as calabash or gourd, is an annual climber belonging to the Cucurbitaceae family. Studies on bottle gourd OVATE genes are limited. In this study, we performed genome-wide identification of the OVATE gene family in bottle gourd, and identified a total of 20 OVATE family genes. The identified genes were unevenly distributed across 11 bottle gourd chromosomes. We also analyzed the gene homology, amino acid sequence conservation, and three-dimensional protein structure (via prediction) of the 20 OVATE family genes. We used RNA-seq data to perform expression analysis, which found 20 OVATE family genes to be differentially expressed based on spatial and temporal characteristics, suggesting that they have varying functions in the growth and development of bottle gourd. In situ hybridization and subcellular localization analysis showed that the expression characteristics of the LsOVATE1 gene, located on chromosome 7 homologous to OVATE, is a candidate gene for affecting the fruit shape of bottle gourd. In addition, RT-qPCR data from bottle gourd roots, stems, leaves, and flowers showed different spatial expression of the LsOVATE1 gene. The ectopic expression of LsOVATE1 in tomato generated a phenotype with a distinct fruit shape and development. Transgenic-positive plants that overexpressed LsOVATE1 had cone-shaped fruit, calyx hypertrophy, petal degeneration, and petal retention after flowering. Our results indicate that LsOVATE1 could serve important roles in bottle gourd development and fruit shape determination, and provide a basis for future research into the function of LsOVATE1.

Identification and Validation of a Core Single-Nucleotide Polymorphism Marker Set for Genetic Diversity Assessment, Fingerprinting Identification, and Core Collection Development in Bottle Gourd.

Germplasm collections are indispensable resources for the mining of important genes and variety improvement. To preserve and utilize germplasm collections in bottle gourd, we identified and validated a highly informative core single-nucleotide polymorphism (SNP) marker set from 1,100 SNPs. This marker set consisted of 22 uniformly distributed core SNPs with abundant polymorphisms, which were established to have strong representativeness and discriminatory power based on analyses of 206 bottle gourd germplasm collections and a multiparent advanced generation inter-cross (MAGIC) population. The core SNP markers were used to assess genetic diversity and population structure, and to fingerprint important accessions, which could provide an optimized procedure for seed authentication. Furthermore, using the core SNP marker set, we developed an accessible core population of 150 accessions that represents 100% of the genetic variation in bottle gourds. This core population will make an important contribution to the preservation and utilization of bottle gourd germplasm collections, cultivar identification, and marker-assisted breeding.

Unraveling the Genetic Architecture of Two Complex, Stomata-Related Drought-Responsive Traits by High-Throughput Physiological Phenotyping and GWAS in Cowpea (Vigna. Unguiculata L. Walp).

Drought is one of the most devasting and frequent abiotic stresses in agriculture. While many morphological, biochemical and physiological indicators are being used to quantify plant drought responses, stomatal control, and hence the transpiration and photosynthesis regulation through it, is of particular importance in marking the plant capacity of balancing stress response and yield. Due to the difficulties in simultaneous, large-scale measurement of stomatal traits such as sensitivity and speed of stomatal closure under progressive soil drought, forward genetic mapping of these important behaviors has long been unavailable. The recent emerging phenomic technologies offer solutions to identify the water relations of whole plant and assay the stomatal regulation in a dynamic process at the population level. Here, we report high-throughput physiological phenotyping of water relations of 106 cowpea accessions under progressive drought stress, which, in combination of genome-wide association study (GWAS), enables genetic mapping of the complex, stomata-related drought responsive traits "critical soil water content" (θcri) and "slope of transpiration rate declining" (KTr). The 106 accessions showed large variations in θcri and KTr, indicating that they had broad spectrum of stomatal control in response to soil water deficit, which may confer them different levels of drought tolerance. Univariate GWAS identified six and fourteen significant SNPs associated with θcri and KTr, respectively. The detected SNPs distributed in nine chromosomes and accounted for 8.7-21% of the phenotypic variation, suggesting that both stomatal sensitivity to soil drought and the speed of stomatal closure to completion were controlled by multiple genes with moderate effects. Multivariate GWAS detected ten more significant SNPs in addition to confirming eight of the twenty SNPs as detected by univariate GWAS. Integrated, a final set of 30 significant SNPs associated with stomatal closure were reported. Taken together, our work, by combining phenomics and genetics, enables forward genetic mapping of the genetic architecture of stomatal traits related to drought tolerance, which not only provides a basis for molecular breeding of drought resistant cultivars of cowpea, but offers a new methodology to explore the genetic determinants of water budgeting in crops under stressful conditions in the phenomics era.

Novel Genomic Regions of Fusarium Wilt Resistance in Bottle Gourd [Lagenaria siceraria (Mol.) Standl.] Discovered in Genome-Wide Association Study.

Fusarium wilt (FW) is a typical soil-borne disease that seriously affects the yield and fruit quality of bottle gourd. Thus, to improve resistance to FW in bottle gourd, the genetic mechanism underlying FW resistance needs to be explored. In this study, we performed a genome-wide association study (GWAS) based on 5,330 single-nucleotide polymorphisms (SNPs) and 89 bottle gourd accessions. The GWAS results revealed a total of 10 SNPs (P ≤ 0.01, -log10 P ≥ 2.0) significantly associated with FW resistance that were detected in at least two environments (2019DI, 2020DI, and the average across the 2 years); these SNPs were located on chromosomes 1, 2, 3, 4, 8, and 9. Linkage disequilibrium (LD) block structure analysis predicted three potential candidate genes for FW resistance. Genes HG_GLEAN_10001030 and HG_GLEAN_10001042 were within the range of the mean LD block of the marker BGReSe_14202; gene HG_GLEAN_10011803 was 280 kb upstream of the marker BGReSe_00818. Real-time quantitative PCR (qRT-PCR) analysis showed that HG_GLEAN_10011803 was significantly up-regulated in FW-infected plants of YD-4, Yin-10, and Hanbi; HG_GLEAN_10001030 and HG_GLEAN_10001042 were specifically up-regulated in FW-infected plants of YD-4. Therefore, gene HG_GLEAN_10011803 is likely the major effect candidate gene for resistance against FW in bottle gourd. This work provides scientific evidence for the exploration of candidate gene and development of functional markers in FW-resistant bottle gourd breeding programs.

Long-read genome assembly and genetic architecture of fruit shape in the bottle gourd.

The bottle gourd (Lagenaria siceraria, Cucurbitaceae) is an important horticultural crop exhibiting tremendous diversity in fruit shape. The genetic architecture of fruit shape variation in this species remains unknown. We assembled a long-read-based, high-quality reference genome (ZAAS_Lsic_2.0) with a contig N50 value over 390-fold greater than the existing reference genomes. We then focused on dissection of fruit shape using a one-step geometric morphometrics-based functional mapping approach. We identified 11 quantitative trait loci (QTLs) responsible for fruit shape (fsQTLs), reconstructed their visible effects and revealed syntenic relationships of bottle gourd fsQTLs with 12 fsQTLs previously reported in cucumber, melon or watermelon. Homologs of several well-known and newly identified fruit shape genes, including SUN, OFP, AP2 and auxin transporters, were comapped with bottle gourd QTLs.

Correction to: A universal pipeline for mobile mRNA detection and insights into heterografting advantages under chilling stress.

[This corrects the article DOI: 10.1038/s41438-019-0236-1.].

A universal pipeline for mobile mRNA detection and insights into heterografting advantages under chilling stress.

Heterografting has long been used to enhance the chilling tolerance of temperature-sensitive crops, including watermelon, whose mechanism is known to involve bidirectional long-distance mRNA movements. Despite several studies reporting on mobile mRNA (mb-mRNA) profiles in plants, accurate identification of mb-mRNAs is challenging owing to an array of technical problems. Here, we developed a bioinformatical pipeline that took most of the known technical concerns into consideration and is considered to be a universal tool for mb-mRNA detection in heterografts. By applying this pipeline to a commercial watermelon-bottle gourd heterografting system, we detected 130 and 1144 mb-mRNAs upwardly and 167 and 1051 mb-mRNAs downwardly transmitted under normal and chilling-stress conditions, respectively. Quantitative real-time PCR indicated a high accuracy rate (88.2%) of mb-mRNA prediction with our pipeline. We further revealed that the mobility of mRNAs was not associated with their abundance. Functional annotation and classification implied that scions may convey the stress signal to the rootstock, subsequently triggering energy metabolism reprogramming and abscisic acid-mediated stress responses by upward movement of effective mRNAs, ultimately leading to enhanced chilling tolerance. This study provides a universal tool for mb-mRNA detection in plant heterografting systems and novel insights into heterografting advantages under chilling stress.

Orphan genes are involved in drought adaptations and ecoclimatic-oriented selections in domesticated cowpea.

Orphan genes (OGs) are genes that are restricted to a single species or a particular taxonomic group. To date, little is known about the functions of OGs in domesticated crops. Here, we report our findings on the relationships between OGs and environmental adaptation in cowpea (Vigna unguiculata). We identified 578 expressed OGs, of which 73.2% were predicted to be non-coding. Transcriptomic analyses revealed a high rate of OGs that were drought inducible in roots when compared with conserved genes. Co-expression analysis further revealed the possible involvement of OGs in stress response pathways. Overexpression of UP12_8740, a drought-inducible OG, conferred enhanced tolerance to osmotic stresses and soil drought. By combining Capture-Seq and fluorescence-based Kompetitive allele-specific PCR (KASP), we efficiently genotyped single nucleotide polymorphisms (SNPs) on OGs across a 223 accession cowpea germplasm collection. Population genomic parameters, including polymorphism information content (PIC), expected heterozygosity (He), nucleotide diversity (π), and Tajima's D statistics, that were calculated based on these SNPs, showed distinct signatures between the grain- and vegetable-type subpopulations of cowpea. This study reinforces the idea that OGs are a valuable resource for identifying new genes related to species-specific environmental adaptations and fosters new insights that artificial selection on OGs might have contributed to balancing the adaptive and agronomic traits in domesticated crops in various ecoclimatic conditions.

Molecular Genetic Mapping of Two Complementary Genes Underpinning Fruit Bitterness in the Bottle Gourd (Lagenaria siceraria [Mol.] Standl.).

Fruit bitterness is a serious problem threatening the bottle gourd (Lagenaria siceraria [Mol.] Standl.) industry worldwide. Previous genetic studies indicated that fruit bitterness in the bottle gourd was controlled by a pair of complementary genes. In this study, based on two non-bitter landraces "Hangzhou Gourd" and "Puxian Gourd," each of which carries a single bitterness gene, and their derived segregation populations, we mapped the complementary genes causing fruit bitterness. Quantitative trait locus (QTL) scanning based on an F2 population detected two QTLs, which was QBt.1 locating in a 17.62-cM interval on linkage group (LG)2 corresponding to a 1.6-Mb region on chromosome 6, and QBt.2 mapped to a 8.44-cM interval on LG9 corresponding to a 1.9-Mb region on chromosome 7. An advanced bulked segregant analysis (A-BSA) well validated the QTL mapping results. Sequence-based comparative analysis showed no syntenic relationship between QBt.1/QBt.2 and the known bitterness genes in cucumber, melon, and watermelon, suggesting that causal genes underlying QBt.1 and QBt.2 were not direct orthologs of the reported cucurbit bitterness genes. Our results shed light on the molecular genetic mechanisms underlying fruit bitterness in the bottle gourd and is useful to guide breeders to properly select parental lines to avoid the occurrence of bitter fruits in breeding programs.

Fine mapping Ruv2, a new rust resistance gene in cowpea (Vigna unguiculata), to a 193-kb region enriched with NBS-type genes.

KEY MESSAGE: A new rust resistance gene Ruv2 was fine-mapped in cowpea to a 193-kb region on chromosome 2, which harboured 23 predicted gene models enriched with NBS-type genes. ZN016 is a landrace vegetable cowpea highly resistant to rust. Two previous studies using mixed-spores inoculation suggested different modes of inheritance of rust resistance in ZN016. In this study, we initially developed a detached leaf assay with a purified single-rust isolate (Auv-LS). Using this approach, we assessed the inheritance of rust resistance in a recombinant inbred line (RIL) population and an F2 population, both derived from the cross of "ZN016" and the susceptible cultivar "Zhijiang282." A single dominant gene mode against Auv-LS was revealed in both populations. QTL mapping showed that this gene was coincident with the Ruv2 locus on LG7, one of the three resistance QTLs previously mapped based on mixed-spores inoculation data. Therefore, Ruv2 was considered as specifically against the rust isolate Auv-LS. Through an analysis of the RIL recombinants at Ruv2, we fine-mapped the gene to an ~ 0.45-cM interval between SNP markers 2_09656 and 2_00973, which corresponded to an ~ 193-kb region on chromosome 2 that harboured 23 predicted gene models enriched with NBS-type genes. Re-sequencing of the two parents revealed polymorphisms in four genes predictively to cause substantial protein structural changes, rendering them valuable candidate genes for future validation. Cross-species syntenic analysis indicated that Ruv2 may represent a novel rust resistance gene in food legumes. A cleaved amplified polymorphic sequences marker tightly linked to Ruv2 was developed to facilitate breeding. This work establishes a basis for map-based cloning of Ruv2 and breeding for rust resistance in cowpea and other legume crops.

GourdBase: a genome-centered multi-omics database for the bottle gourd (Lagenaria siceraria), an economically important cucurbit crop.

GourdBase is an integrative data platform for the bottle gourd to examine its multifarious intuitive morphology and annotated genome. GourdBase consists of six main modules that store and interlink multi-omic data: the genome (with transcriptomic data integrated) module, the phenome module, the markers/QTLs module, the maps (genetic, physical and comparative) module, the cultivars module, and the publications module. These modules provide access to various type of data including the annotated reference genome sequence, gene models, transcriptomic data from various tissues, physical and comparative genome maps, molecular markers in different types, phenotypic data for featuring traits including fruit shape and umami taste, and quantitative trait loci (QTLs) that underlie these traits. GourdBase is intuitive, user-friendly and interlinked and is designed to allow researchers, breeders and trained farmers to browse, search and fetch information on interests and assist in genomics-driven studies and breeding. The knowledge base and web interface can be accessed at http://www.gourdbase.cn/ .

Genomic regions, cellular components and gene regulatory basis underlying pod length variations in cowpea (V. unguiculata L. Walp).

Cowpea (V. unguiculata L. Walp) is a climate resilient legume crop important for food security. Cultivated cowpea (V. unguiculata L) generally comprises the bushy, short-podded grain cowpea dominant in Africa and the climbing, long-podded vegetable cowpea popular in Asia. How selection has contributed to the diversification of the two types of cowpea remains largely unknown. In the current study, a novel genotyping assay for over 50 000 SNPs was employed to delineate genomic regions governing pod length. Major, minor and epistatic QTLs were identified through QTL mapping. Seventy-two SNPs associated with pod length were detected by genome-wide association studies (GWAS). Population stratification analysis revealed subdivision among a cowpea germplasm collection consisting of 299 accessions, which is consistent with pod length groups. Genomic scan for selective signals suggested that domestication of vegetable cowpea was accompanied by selection of multiple traits including pod length, while the further improvement process was featured by selection of pod length primarily. Pod growth kinetics assay demonstrated that more durable cell proliferation rather than cell elongation or enlargement was the main reason for longer pods. Transcriptomic analysis suggested the involvement of sugar, gibberellin and nutritional signalling in regulation of pod length. This study establishes the basis for map-based cloning of pod length genes in cowpea and for marker-assisted selection of this trait in breeding programmes.

Natural variation and gene regulatory basis for the responses of asparagus beans to soil drought.

Asparagus bean (Vigna unguiculata ssp. sesquipedalis) is the Asian subspecies of cowpea, a drought-resistant legume crop native to Africa. In order to explore the genetic variation of drought responses in asparagus bean, we conducted multi-year phenotyping of drought resistance traits across the Chinese asparagus bean mini-core. The phenotypic distribution indicated that the ssp. sesquipedalis subgene pool has maintained high natural variation in drought responses despite known domestic bottleneck. Thirty-nine SNP loci were found to show an association with drought resistance via a genome-wide association study (GWAS). Whole-plant water relations were compared among four genotypes by lysimetric assay. Apparent genotypic differences in transpiration patterns and the critical soil water threshold in relation to dehydration avoidance were observed, indicating a delicate adaptive mechanism for each genotype to its own climate. Microarray gene expression analyses revealed that known drought resistance pathways such as the ABA and phosphate lipid signaling pathways are conserved between different genotypes, while differential regulation of certain aquaporin genes and hormonal genes may be important for the genotypic differences. Our results suggest that divergent sensitivity to soil water content is an important mechanism configuring the genotypic specific responses to water deficit. The SNP markers identified provide useful resources for marker-assisted breeding.

Association Mapping for Fusarium Wilt Resistance in Chinese Asparagus Bean Germplasm.

Fusarium wilt (FW) is an important vascular disease attacking asparagus bean [Vigna unguiculata (L.) Walp. subsp. unguiculata Sesquipedalis Group] in China. The level and genetic variability of FW resistance in the Chinese asparagus bean germplasm remains elusive. In the current study, FW resistance was assessed across a natural population consisting of 95 asparagus bean and four African cowpea [Vigna unguiculata (L.) Walp. subsp. unguiculata Unguiculata Group] accessions. The disease index (DI) based on the severity of leaf damage (LFD) and vascular discoloration (VD) varied highly across the population and the highly resistant varieties used for vegetable are very limited. Genome-wide association study identified 11 and seven single nucleotide polymorphisms (SNPs) that are associated with LFD and VD traits, respectively. These SNPs were distributed on nine linkage groups of the asparagus bean genome and each accounted for less than 5% of the phenotypic variation. Overall, the nonstandard vegetable (NSV) subgene pool harbors favorable alleles in a higher frequency than the standard vegetable (SV) subgene pool. Individual NSV-type accessions tend to possess a greater number of favorable alleles than the SV-type ones. A SNP marker 1_0981 was converted to a cleaved amplified polymorphic sequences (CAPS) marker to facilitate future breeding. To our knowledge, this is the first report of an association mapping (AM) study in asparagus bean. The results obtained suggests that resources for FW resistance is relatively limited in the SV subgene pool; hence, introducing resistant alleles from the NSV accessions into currently leading SV cultivars will be important to improve FW resistance of the latter.